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We compute quantum and thermal corrections to the right-handed 
neutrino interaction rate in the early universe at next-to-leading order 
in all the relevant SM couplings (gauge, top Yukawa and higgs cou- 
plings). Previous computations considered 2 — > 2 scatterings, finding 
Ph infra-red divergences. The KLN theorem demands that infra-red di- 

Pu vergences cancel in the full result: after adding 1 — > 3 and one-loop 

^ virtual corrections that enter at the same order we find a simple result. 
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1 Introduction 



Thermal leptogenesis [1] seems the most plausible explanation of the observed baryon asymme- 
try of the universe [2] . A key quantity for leptogenesis is the (space-time density of the) rate at 
which the thermal plasma of the early universe at temperature T creates quanta of the lightest 
right-handed neutrino N with mass M: 

( rp\ dN N 

7 " (T) = dVdt (1) 

In thermal equilibrium, the creation rate equals the destruction rate, such that both quantities 
are usually named "equilibrium interaction rate" . It enters in the Boltzmann equation for the 
evolution of the total N abundance n^'- 

dY N (Y N \ 
sHz —r- = ~ Utsq ~ 1 7jv, (2) 



dz ; 

where Y/y = n^/s, s is the entropy density, z = M/T and H(z) is the expansion rate. 

At leading order, 7jv is given by the thermal average of the N — » LH, LH decay rate T^{E) 
induced at tree level by the Yukawa coupling A NLH, where L and H are the usual lepton and 
Higgs doublets and E = ^Jp 2 + M 2 is the energy of N [3, 4, 5, 6]: 



^ = 2 J (#^ r - Tn {e) = A 8^¥' (3) 

where /n = l/(e E ^ T + 1) is the Fermi-Dirac distribution. 

The goal of this paper is computing all quantum and thermal corrections to the N interaction 
rate 77V, up to NLO in all relevant SM couplings g: the gauge couplings gi,gy, the top Yukawa 
coupling X t = m t /v and the Higgs quartic interaction X h = (m h /v) 2 , where v = 174 GeV. 

Previous partial results are extremely complicated because only some NLO effects have been 
computed, missing the great simplification that happens when including all NLO corrections: 
infra-red (IR) divergences cancel out in the total result, as demanded by the Kinoshita-Lee- 
Nauenberg (KLN) theorem [7]. The final result must have the form 



77V = 7^° 



1 + + K T g^ 2 + 0{^f + NNLO orders] , (4) 

where K Q and Kt are order-one constants, computed in the rest of this paper. 

• Kq is the zero temperature quantum correction, that was so far ignored. 

• Kt is the finite temperature correction, that is precisely needed only at T M. Indeed 
thermal leptogenesis does not depend on the initial iV abundance only if A is large enough 
that right-handed neutrinos remain close to thermal equilibrium down to low temperatures 
T M; we therefore only need a precision computation of 7jv in such limit. 

Some thermal effects have been computed so far [4, 5, 8, 9, 6]: the contribution coming from 
2—7-2 scatterings (such as AN — > LH, where A is any SM vector), finding lengthy expressions 
where thermal masses regulate infra-red (IR) divergences [8, 5]. In thermal field theory this 
effect is just one correction to 7^; after adding all other effects (3-body decays, such as iV — > 
LHA, and virtual corrections) we will find that IR divergences cancel out. 

In section 2 we compute quantum corrections. In section 3 we compute thermal corrections. 
The computations are lengthy, but thanks to cancellation of IR divergences the final result can 
be written in one line: it is presented in the conclusions, section 4. 
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2 Quantum corrections 



In this section we compute the quantum correction to the N interaction rate T^, up to 
0{g\, <7y, , Aft). Such corrections have been neglected so far, and include two effects: i) one 
loop corrections to N — » LH and ii) 3-body decays, such as iV — )■ LHA. Separately they are 
infra-red divergent. However, as computed in the rest of this section and as demanded by the 
KLN theorem [7], they combine to produce a N decay rate that does not depend on the IR 
structure of the theory. 



2.1 Quantum corrections: tools 

We employ dimensional regularization for both IR and UV divergences. The phase space in 
d = 4 — 2e dimensions for one particle with quadri-momentum P = (E,p) decaying into n 
particles with quadri-momenta Pj is: 

d Pi n rr j-,2 2\ d Pi 

upi, api = ' ■ - 

i=l i=l 

For massless final-state particles the 2-body phase space is 



d$ n = {2^5{p - y: Pi) n ®, ^ - 6 ^ - ^ = (2^%- (5) 



M rf - 4 r(d/2 - 1) ^4 }_ 
2 2 d ~ 1 n d / 2 - 1 T(d - 2) 8tt { } 

and the 3-body phase space is 

rf$3 = 16(2vr)3 ^ J f{d^2) dXl dX2 = 128.3 ' (?) 

where = 2P, • P/P 2 such that xi + x 2 + x% = 2. The integration region is < x\ < 1 and 
1 — xi < x% < 1: it is obtained considering the triangle with sides xi,X2,X3 and demanding 
that any side is longer that the difference of the other two and shorter than their sum. 

We now compute the relevant corrections, in increasing order of difficulty: higgs, top and 
gauge. 

2.2 Higgs quantum correction 

There are no corrections induced at one loop by the quartic higgs coupling A/j. 



2.3 Top quantum correction 

We compute the quantum corrections induced by the top quark Yukawa coupling X t HQU. 
The only virtual correction is the correction to the H propagator. All particles in the loop 
are massless, such that this correction vanishes in dimensional regularization. The only NLO 
correction is then the 3-body decay iV — > LQU . We find: 

\2 / o 91 \ 

T{N^IS)V)=T^{-- + U--y (8) 

where here and in the following i = InM 2 //* 2 , with /2 being the MS renormalization scale. 
The UV divergence gets reabsorbed by writing the top Yukawa coupling in terms of its value 
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renormalized in the MS scheme and the scale (i = M, obtaining, for the top quantum corrections 
at NLO: 



, top>T=0 _ A 2 (M)M 



A' 



S7T 



1 



21 Jl 
2 (4?r) 



(9) 



Since IR divergences cancel, the same result can be obtained with different IR regularizations. 1 
2.4 Gauge quantum corrections 

We consider one abelian vector with coupling a under which L and H have charge one; it will 
be easy to add at the end the group factors appropriate for SM vectors. Virtual corrections to 
on-shell propagators of massless particles vanish in dimensional regular ization. Only the vertex 
diagram contributes to virtual corrections, and the result is: 



virtual 



(N ->• LH) = r 



a. 



0~ 



4 I - 1 77T 2 



(12) 



A-n \ e 2 e v '3 

Emission of one vector A from either the fermion L or the scalar H gives, in Feynman gauge: 

a 



T(N LEA) = T 



47r \e 



4 7-4£ a . n , 7tt 2 45 

+ + - 7) - — 
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(13) 



Summing real and virtual corrections the IR divergence cancels, leaving 



NLO 



To 



29 



a ( 3 u 2 



MA 2 (M) 



8tt 



1 + 



29 a 
T 4^ 



(14) 



having expressed T in terms of A renormalized at \x — M in the MS scheme. The result is IR 
convergent as guaranteed by the KLN theorem. Finally inserting the group factors appropriate 
for the SM electroweak vectors, a — > (ay + 3a 2 )/4, we find: 



r gauge,T= 
1 iV 



=0 MA 2 (M) 



S7T 



1 + 



29 
327 



(3a 2 + «y) 



(15) 



Since IR divergences cancel, the same result can be obtained with different IR regularizations. 



1 Alternatively, one can regularize the IR divergence with a small top mass, m t < mjv- In such a case we 



get 



r(jv -> lqu) = r 



23 m t 

■ — — 6 m 
2 M 



The virtual correction is both UV and IR divergent: 



(10) 



rtual 



(N -► Li/) = r 



(4tt) : 



(11) 



giving again the same final result for T^ P ' T = r + T(N — > LQU) + r virt uai(^V — > LH). 
2 Using a small vector mass vtia as infrared regulator we find 



IYTV LHA) = T x - 2vr 2 +421nr + 241nr 2 ), r = — < 1. 

127T 2 rriTv 



(16) 



The correction to the kinetic terms of a massless fermion and of a massless scalar due to a loop of an abelian 
vector with mass my in a generic £ gauge are 



PPl 



a ,.,\ u 2 . , 3, 
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2.5 Summary 



Including quantum corrections at NLO in all relevant SM couplings, the N decay rate does not 
receive any IR-divergent correction and is: 



r T=0 
1 N 



MX 2 (M) 



8tt 



29 ,o \ 21 A 2 



32tt 



2 (4tt) 



(17) 



where A(M) is the neutrino Yukawa coupling renormalized at the MS scale \i = M. It satisfies 
the well known RGE equation: 



(4tt 



dlnfl 



-A 



;9 Y 



at - 3A, 2 



RGE equations at LO have been computed in [10], but the connection between M and A(M) 
with neutrino masses has not yet been computed at NLO. 
In the next section we consider thermal corrections. 



3 Thermal corrections 

As already discussed in the introduction, we are interested in the dominant thermal corrections 
in the low temperature limit, T <C M. 

At tree level, thermal corrections are exponentially suppressed by the Boltzmann factor 
e~ M ' T , as clear by the thermal function in (3). Such corrections have been included in 
previous works [5]. 

At loop level, thermal corrections are only power suppressed (e.g. because in processes such 
as iV — > LHA the vector energy can be comparable to the temperature even at M ^> T): we 
want to compute the dominant corrections proportional to (T/M) 2 . Despite the (T/M) 2 sup- 
pression, thermal corrections are relevant because not suppressed by any l/(47r) 2 loop factor, 
as anticipated in eq. (4), where the coefficient is expected to be of order one. Inserting 
such interaction rate into Boltzmann equations, we find the correction to the baryon asym- 
metry shown in fig. 1, equal to — 1.5%K T (—3.5%K T ) for rh\ = m atm (m sun ), so that thermal 
corrections are expected to be comparable to quantum corrections and need to be computed. 

We compare our approach with previous works. 

• Ref. [5] performed a resummation of the leading thermal corrections in the high-temperature 
regime, T > M, approximatively described by thermal masses: 

m 2 H = (^2 + ^y + ^ 2 + ^A^T 2 , (18a) 

^ .2 , 1 2 \ th2 



< = (32* + 32ftj T - < 18b > 

We are interested in the low-temperature regime, that is the relevant one in the regime 
rhi ^> 10~ 3 eV where leptogenesis is computable independently of the initial iV abundance. 



Adding also the vertex diagram, the total virtual correction is ^-independent (such that RGE equations for the 
couplings are gauge-independent): 



virtual 



°L ( 3 + 3 jn _ Minr - 81n 2 r + — 
4tt I e M 2 3 



5 




rh\ in eV 

Figure 1: Correction to the baryon asymmetry due to thermal effects at tree level (red line, 
suppressed by e~ M ' T ) and at loop level (blue line, suppressed by T 2 /M 2 ). The shaded region 
shows the uncertainty due to the initial right-handed neutrino abundance, varied from negligible 
to dominant. The dashed lines show the values rh\ = (Am^) 1 '' 2 and (Am 2 tm ) 1 '' 2 . 



The thermal masses 3 m 2 H and 2m 2 L of L and H reduce the phase space for iV — > LH decays, 
providing one contribution to Kt = —2m 2 H /T 2 = — 3(7|/8 H ~ —0.2. 

• Another contribution to comes from 2 — > 2 scatterings, such as AN — > LH. This 
process reduces to iV — > LH in the limit of vanishing energy of the vector A, and it gives 
an IR divergent contribution to Kt (it was computed in [5] using thermal masses as IR 
regulator). Indeed such scattering rate is proportional to the number density of the initial 
vectors, given by the Bose-Einstein distribution ua(uj) = 1/[1 — e w / T ], that diverges at 
small vector energy u — > 0, giving rise to the IR divergence. 

These results are IR-divergent because incomplete: e.g. IR-divergent terms proportional to ua 
cancel after adding AN — > LH scatterings with N — > LHA decays and with virtual corrections 
to iV — >• LH. 



3.1 Thermal corrections: tools 

We want to compute thermal corrections induced by the large gy,92, ^t, couplings, while we 
can neglect those induced by the smaller neutrino Yukawa interaction A. In such a case the 
interaction rate at finite temperature is precisely defined and computed from the imaginary 
part of the iV propagator in the thermal plasma, computed by explicitly summing all possible 
cuttings of the relevant Feynman diagrams shown in fig. 2 and fig. 5: the first one describes 
the tree-level result. 

3 The factor of 2 in front of m\ arises because thermal masses are conventionally defined at zero momentum 
p (energy of a particle at rest with respect to the plasma) . Thermal dispersion relations are not relativistic, and 
for fermions the "thermal" mass squared at p 3> T is 2m 2 L . 
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According to the real-time formalism of thermal field theory [11], the decay/absorption rate 
T of a particle with mass M and quadri-momentum P = (E,p, 0,0) coupled to a thermal 
plasma by a weak coupling A is given by, at leading order in A 

T(E> = (19) 

where II is the propagator (if the particle is a boson), or its spin-average (if the particle, as 
in our case, is a fermion: H(P) = Tr[(.p + M)£(P)]/2 where £ is the fermion propagator). 
Il > is the non time-ordered propagator, and can be computed following the rules of Kobes 
and Semenoff for imaginary parts of Feynman diagrams at finite temperature, which generalize 
the cutting rules valid at T = for imaginary parts of Feynman diagrams. II > is essentially 
equivalent to Imll and gives rise to simpler expressions. 

The two-point function U > (xi,X2) is computed summing over all possible ways of circling 
the internal type I Zi vertices; the external vertex X\ is circled and the external vertex x 2 is 
uncircled. Non time-ordered correlations, such as II > , are computed using Feynman diagrams 
with the following additional rules: 

• Reverse the sign of the coupling associated to a vertex, if it is circled. 

• For a propagator connecting two points x and y: 

i) Use the standard propagator P(x — y) if neither x nor y are circled; 

ii) Use the propagator P* (x — y) if both x and y are circled; 

iii) Use the propagator P < {x — y) if x but not y is circled; 

iv) Use the propagator P > (x — y) if y but not x is circled; 

In momentum space, the propagators P = A B for scalars, P = (Ifc + m)Ap for fermions (when 
P* is needed, the complex conjugate does not act on 7^ matrices), P = —g^ v ls. B for vectors in 
the Feynman gauge are 

A B = A + 27m B (K )5(K 2 - m 2 ), A F = A - 27m F (K )5(K 2 - m 2 ), 

A> = [9(+K ) + n B (K )]2n5(K 2 - m 2 ), A> = [9(+K ) - n F (K )}2Ti5(K 2 - m 2 ), (20) 

A< = [9{-K ) + n B (K )]2n5(K 2 - m 2 ), A< = [9{-K ) - n F (K )}2n5(K 2 - m 2 ), 

where K is the quadri-momentum and A = i/{K 2 — m 2 + ie) is the propagator at T = 0; A > 
(A < ) are applied when K enters into (exits from) the circled vertex. The functions n BtF (E) = 
f BF (\E\) are the Bose-Einstein and Fermi-Dirac statistical factors, f B (E) = l/[e E ^ T — 1] and 
f F (E) = l/[e E / T + l}. 

After lengthy manipulations one finds that ill-defined products of 5 functions cancel out 
when summing over circlings of each type of diagram separately, and one recovers a decompo- 
sition into 'real' and 'imaginary' contributions and expressions similar to the ones well-known 
at T = 0. 

Thermal corrections can be written in a way similar to quantum corrections at T = by 
replacing propagators with thermal propagators and phase space with 'thermal phase space'. 
For a particle with mass M and quadri-momentum P this is defined as: 

dP = dp [9(P ) ± n B>F (P )} = ~- d [9(P ) ± n BjF (P )} 2n5(P 2 - M 2 ), (21) 
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which generalizes the T = phase-space integral in d = 4 — 2e dimensions. Thermal processes 
such as the N interaction rate combine contributions from different quantum processes (e.g. 
decays such as iV — > LHA and scatterings such as AN — » LH): they are unified by the 'thermal 
phase space' in (21) that allows each particle to have positive Pq (which means it is in the final 
state, as usual) and negative Pq (which means it is in the initial state). Indeed the expression 
in square brackets gives the statistical factors appropriate for the two cases: 1 ± Ub,f i n the 
final state and in^^ in the initial state. At T = the thermal factors ns t F vanish and the 
^-function forces Po = E > 0, recovering the usual phase space for final state particles: 

lim dP = dp= , ^ f - . (22) 
t^o F {2n) d ~ l 2E y ' 

For example, at leading order only the processes N — > LH, LH are present 

^ (E) = ^= fdP L dP H ■ (2n) d 6(P -P H -P L ).^ Mf , (23) 

J final 

where the sum is over polarizations and quantum numbers in all final states. Kinematics 
demands that L and H can only be in the final state, such that the thermal average of such 
decay rate is the usual expression [5]: 

7 £» = j dp N dp L dp H f N (l + f H )(l - f L )(2n) d 5(P -P L - P H ) (24) 

J all 

where the sum is over initial and final polarizations and quantum numbers. Thermal effects 1±/ 
break the usual Lorentz dilatation of the decay rate, T(E) = T(M)M/ E, such that the integrals 
cannot be all done analytically. However these thermal effects are irrelevant at T <C M, because 
suppressed by Boltzmann e -°( M )/ T factors, and the tree result can be approximated in terms 
of the decay width of a iV at rest at T = 0: 

7 r^ M = M)xH^> (25) 

where the ratio of standard Bessel functions Ki/K 2 — 1 is the thermal average of T(E) /T{M), 
approximated as M/E according to the standard relativistic formula. 

At NLO order in the gauge couplings one has an extra vector A with quadri-momentum K 
that enters in various 3 — » 1 and 2 — > 2 processes, all described by the following expression 
that resembles the T = result: 

IT 1 = ^ J dP L dP H dK ■ \JZ\ 2 ■ (2n) d 5(P -P l -Ph- K), (26) 

where the phase-space factors dP include the statistical factors for each particle, see eq. (21). 
The L, H, A particles can be in the initial or in the final state, according to the sign of 
P°, Pff, K°: thereby eq. (26) unifies the 1 — > 3 decay with 2—^-2 scatterings. 

We see that we only need to compute at NLO the decay width of a N at rest with 
respect to the thermal plasma. 

When computing thermal corrections we will also get (in the T — > limit) the quantum 
corrections, already computed in section 2. They must be discarded keeping only purely thermal 
corrections, which are not affected by UV divergences. 
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Figure 3: Imaginary parts of the relevant Feynman diagrams for the top corrections, using 
the Kobes and Semenoff circling notation [11]. When circled vertices form connected sets, we 
also show the corresponding cutting notation that holds at T = 0. Momenta indicated in the 
diagrams flow from left to right. 



3.2 Higgs thermal correction 

The only thermal effect present at NLO is the higgs coupling contribution to the higgs thermal 
mass, see eq. (18a). It reduces the decay width for iV — >■ LH by a factor 1 — 2m 2 H /M 2 , such 
that the higgs coefficient in eq. (4) is 

^ iggs = -1. (27) 



3.3 Top thermal correction 

We compute here the thermal corrections induced by the top quark Yukawa coupling A t HQU 
at 0(T 2 /M 2 ). We will find that the top coefficient in eq. (4) is 

= 0. (28) 

Some of the formulae introduced in this section will be used to compute the gauge corrections 
in section 3.4. 
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3.3.1 Virtual top thermal correction 

Let us first consider the virtual top corrections, diagrams el, e2 and e4 in fig. 3. The last 
diagram is proportional to an ill-defined 5 function squared, that cancels exactly with the 5 
function squared coming from the other diagrams (including the real correction e3). We are 
therefore left with the remaining part of diagrams el and e2: thermal self energy corrections 
to the H line. Their contribution can be written as: 

ZT - 1 - , (29) 

where m^ top and Z l ^ p are the top Yukawa contribution to the thermal mass and thermal 
wave function renormalization constant of H. These quantities can be computed from the top 
contribution to the thermal H self-energy: 

U t (P H ) = h\W~ d Tr J ^ d S( P Q) S ( P Q ~ P ")> (30) 

where S(Pq) is the thermal propagator of the quark, of which, as before, we neglect the mass. 
What actually appears in the calculation is not the full H self-energy, but only its real part, 
Re IT. Since Lorentz invariance is broken at finite temperature, Il f does not only depend on Pjj 
but separately on Pfj and pfj. The thermal mass however turns out to be Lorentz invariant, 
and can be computed from ReIL(0), reproducing the standard result of eq. (18a): 

<to P = 6A? f ^{P 2 Q )n F {P m ) = \%T\ (31) 
The wave function renormalization is given by the next to leading term in the expansion around 

ZT = l + ?^{f H , p 2 ), (32) 

which generalizes the T = formula. We obtain = 1. So the virtual top correction to the 
higgs decay width is just due to the reduction in phase space arising from the thermal Higgs 
mass [5]. For T<CM we have: 

rS-^C^O = r tree (£) (-yw) le - ^ op ' virtual = ~\ ^ 

3.3.2 Real top thermal corrections 

The real correction (diagram e3 in fig. 3) leads to the decay /absorption rate as described in 
section 3.1 for the process iV — >■ LQU, LQU with quadri-momenta P = Pl + Pq + Pu- 

rlZ(E) = ^ J dP L dP Q dP v ■ |^ top | 2 • (2n) d 5(P - P L - P Q - P v ), (34) 

where 

l^topl 2 = 3A 2 A?^^. (35) 
Pq ■ Pu 

The 3-body phase space can be simplified as follows. We integrate over dPjj using conservation 
of energy-momentum. Next: 

5{P V ) = 5((P -P L - P q ) 2 ) = 5(P 2 - 2P ■ (P L + P Q ) + 2P L ■ P Q ). (36) 
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From now on, as previously discussed, we can assume that P is at rest with respect to the 
plasma, P = (M, 0,0,0); then the statistical functions n F depend on Plo, Pqo and Puo only 
and (36) becomes 



5{M 2 - 2M(P L0 + Pqo) + 2P L0 Pq - 2\P L0 P Q0 \ cos 9), (37) 

where Plo and Pqo can be positive or negative and 9 is the angle between j>l and pq. The 
phase-space integral over 9 is done using the 5 function of eq. (37). The condition | cos6*| < 1 
gives the allowed regions, that we now explicitly compute in terms of the two relevant free 
variables. 

Common choices are P^o and Pqo, directly related to the dimensionless variables already 
employed in the quantum computation: xl,u,q = 2Pl,u,q 4 P '/ P 2 = 2Pl,u,qo/M (xl + xu + xq = 
2). It is more convenient to use two slighly different variables z and y, defined by: 

x Q = z, xu = 1 - yz, x L = 1 - z(l - y) (38) 

16 M M, M, 

Pqo = y z > P/o = y(l-^), Plo= (39) 

Inverting these relations: 

, = XQ ^2^, , = i + _=i + _ =! + ___. (40) 

The z,y variables determine the angle as | cos^l = |1 — 2y/[l — (1 — y)z]\ such that the con- 
dition | cos 6*| < 1 gives the kinematically allowed regions, plotted in fig. 4a in terms of the 
dimensionless variables z, y. We just rediscover the well known allowed 1 — > 3 decay and 2 — > 2 
scatterings. 

We can now rewrite the 3-body phase space of (7), in terms of the z,y variables: 

d ® 3 ~ 16(2vr)3 ) f{d~2) dz dy ~ 128*3 " (41) 

Also the amplitude of eq. (35) can be rewritten in terms of y and z: 

\^ top \ 2 = 3X 2 X 2 t 1+ {1 iy ^ Z . (42) 

By considering separately the various processes corresponding to the four integration regions in 
the left plot of fig. 4, we finally obtain the real top corrections as the sum of their contributions: 

T n ^lqu = r t ° ree -r^ / dy dz(l-n F (PLo))(l-n F (P QO ))(l-n F (P uo ))F(y,z),(43&) 

o \2 poo /*oo 

T nu ^lq = -rieej^l dyj dz(l-n F (PLo))(l-n F (P Q o))n F (P uo )F(y,z), (43b) 

Q\2 pi rO 

Tnq^lu = -r t ° rcc — ^ / dy dz(l-n F (PLo))n F (P Q o)(l-n F (P uo ))F(y,z), (43c) 

i07r JO J-oc 
q\2 rO poo 

Tnl^qu = -r? ree — ^ / dy dzn F (P LO )(l-n F (P Q0 ))(l-n F (P U0 ))F(y,z), (43d) 

1D7T J_ OQ J x 
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top corrections 



2 



NL->QU 



L collinear to U 



N^LQU 
soft Q 



NQ->LU 



NU^LQ 



1 

y = i + (P w -M/2)/K 



gauge corrections 



NL->HA 



L collinear to // 



- 



N->LHA 

soft A 



NA^LH 



NH^LA 



1 

v = 1 + (P w -M/2)/K 



(a) (b) 

Figure 4: Kinematically allowed regions and their boundaries, where one has soft and collinear 
singularities (thick lines) for top Yukawa (left) and gauge (right) thermal corrections. Only the 
decay process is present at T = 0. 



where r° ree is the tree level decay rate at T = and 



2 \d/2-2 



r(d/2-i) (i-y)^ 



We can now check that infrared divergences cancel after summing the different processes. The 
situation is relatively simple: 

• There are no thermal IR singularities, because we only have fermionic distributions, which 
are not singular at n F (0); 

• There is no soft singularity at z — > 0: the infrared singularity in |^top| 2 gets canceled by 
positive powers of z in the 3-body phase space in (41); 

• There is a collinear divergence, described by the 1/(1 — y) factor at y — > 1. As clear 
from fig. 4 this corresponds to Q collinear with U and the first 3 processes in eq. (43) are 
separately infra-red divergent. 

One infrared divergent process is the decay iV — > LQU which was not included in previous 
works, that thereby missed the cancellation of infrared divergences, which takes place in the 
full result, as we now describe. 
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The quantities in (43a)-(43c) have the following poles at d — 4: 



j; 1 MX 2 3A? r ,M X1 Z" 1 , r ,Mz„ ,M(l-z) N1 , 

„div 1 A/A 2 3A 2 ,M,,r,„ ,Mz., ,M(l-z), ,, r , . 

rtiv 1 MA 2 3A| . M f" , ,Mz. u ,M(l-z).. 

The integrations over z can be performed analytically, with the result: 

r#W + r^zQ + ^ = ^r^^^t 1 - My)][i + My)]- (46) 

In the limit T — > (i.e. neglecting thermal functions) we recover the infra-red divergent real 
quantum correction (already computed in section 2.3), here contained in F^Z+lqu- For generic 
T, integrals over fermionic thermal functions squared produce the bosonic thermal function 
1 + tib(M/2), such that the total real infra-red divergence in real processes has the same 
structure of the tree level result. 

Such real IR divergence cancels out with the IR divergence in the virtual contribution. The 
virtual contribution vanishes because it has opposite IR and UV divergences, not distinguished 
by dimensional regularization. So, formally, the cancellation takes place when expressing the 
tree-level result in terms of the renormalized coupling: 

A2 = ^ ( 1 " + -) ■ < 47 > 

where A(M) is, as before, the MS coupling at the scale JJl = M and the dots are the contributions 
of the other relevant coupling constants. Let us see in detail how this occurs. The tree level 
contribution at finite temperature can be written as 

r l „.(M) = — (^j r (d -2) V MylHl+My)]. (48) 

Substituting (47) into (48) and setting Jx = M one obtains a divergent part which is exactly 
opposite to (46). 

Having verified that IR divergences cancel, the practical recipe to get the finite result is to 
sum the finite parts of the various contributions: 

^N^LQU - FnU-LQU + ^NU^LQ ~ ^NU-^LQ + ^ NQ^LU ~ ^NQ-yLU + ^NL^QU- (49) 

The integrals in (43a)- (43d) can be computed numerically for arbitrary values of T/M and 
analytically for small T/M, following the appendix of [12]. 

We can neglect all terms exponentially suppressed by factors ~ e~ M l T . Each np(E) in 
(43a)-(43d) provide such suppression factor, unless kinematics allows integration regions where 
E <T. Kinematics allows at most one particle among Q,U,L to have energy E < T; thereby 
we can keep only the terms of (43) which have one thermal function (dropping the T = term 
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with no thermal functions and the terms with two thermal functions): 

T N _ LQU ~ _r t ° rce — ^ / dy dz(n F (P L0 ) + n F (P Q0 )+n F (P U0 ))F(y,z), (50a) 
Ion" Jo J 

2^2 /"oo poo 

T NU ^ LQ ~ - r t°ree^ / ^/ dz n F {P U0 )F{y , z) , (50b) 
3A 2 f 1 Z" 

J- NQ^LU — -T° tiee --\ / / ^n F (P Q0 )F(y,z), (50c) 
Tnl^qu ^ -T° tiee --\ dy dzn F (P L0 )F(y,z). (50d) 

1D7T J_ OQ J 1 

Finally, we integrate along the phase space where some particle has E < T: 

• Region Pqq ~ 0, corresponding to the line z ~ (see fig. 4a) The relevant terms are the 
second one of (50a) plus (50c). We expand F(y,z) in powers of z. To this end notice 

(1 - z) d / 2 - 2 [l + {y- l)z] = 1 + (1 - d/2 + y)z + 0(z 2 ). (51) 

The O(z) part on the right hand side of (51) gives no contribution to Kj<: indeed its 
contribution to the second term in (50a) is opposite to that to (50c) up to higher orders 
in (T 2 /M 2 ). So the only contribution at first order in T 2 /M 2 comes from the O(z) part 
and we obtain 

Ar(|p QO |~o)^r° ree |^. (52) 

• Region Puo ~ 0, corresponding to the point y ~ z ~ 1 (see fig. 4a). The relevant terms 
are the third one of (50a) plus (50b). Let us first expand F(y, z) in powers of (z — 1) by 
noticing 

1 ? [1 + (y - l)z] = y + (dy-l-3y)(z-l) + 0{{l - zf). (53) 



The contribution coming from the 0(1) part in (53) vanishes up to 0(T 2 /M 2 ), while the 
Oil — z) part gives 

Ar^ol-o^r^!^. (54) 

• Region P L ~ 0, corresponding to the point y ~ and z ~ 1 (see fig. 4a). The relevant 
terms are the first one of (50a) plus (50d). We find none of these terms contribute up to 
0(T 2 /M 2 ). 

If we now sum (52) with (54) we obtain x^ v ' real = 1/2. So the real plus virtual corrections, eq. 
(33), give the final result anticipated in (28): K^ p = 0. 

3.4 Gauge thermal correction 

We now turn to the thermal gauge corrections to the interaction rate of iV up to 0(T 2 /M 2 ) 
or, in other words, the gauge corrections to K T in eq. (4). Our result is 

Kf nge = 0. (55) 
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Figure 6: Imaginary parts of the relevant Feynman diagrams for the gauge corrections. We use 
the same notation as in fig. 3. 
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3.4.1 Virtual gauge thermal corrections 

Let us start with the virtual corrections. We divide the contributions in vertex correction 
(diagrams dl and d2 in Fig. 6), corrections to the L propagator (diagrams bl, b2 and b4) and 
corrections to the H propagator (cl, c2 and c4). These diagrams separately contain ill-defined 
squares of 5 functions, which, like in the top case, cancel when one sums them with the real 
corrections (diagrams b3, c3, d3 and d4). 

We first consider the virtual contribution from the vertex diagrams dl and d2. Like for 
the real part (see section 3.1), the KS rules again give an expression that resembles the T = 
result: 



vertex 



-^X 2 x^~ 2d J dP L dP H (27r) d 5(P - P L - P H ) x g 2 J ^-JT 9, (56) 

where g 2 = 3gf/4 + gy/4; the first integral is over the 'thermal phase space' of the two-body 
decay; the second integral is the loop integral over the vector quadri-momentum K; JV is the 
usual numerator coming from the Feynman diagram: 

JT = (2P H - K) p Tr[f>f> L y p (f> L + $)} , (57) 

and 8? is the product of the various thermal propagators A BF , that can be expanded in terms 
of zero temperature propagators A as follows: 

& = Re[A F (P L + K)A B (P H -K)A B (K)} = 
= Re[A (P L + K)A (P H -K)A (K)] + 

+2n5(K 2 )n B (K ) Re[A (P L + K) A (P H - K)} + (58) 
+27r5((P H - Kf)n B {P Q H - K°) Re[A (P L + K)A (K)) + 
-2nS((P L + K) 2 )n F (P° L + K°) Re[A (P H - K)A (K)]. 

All other terms with higher powers of n B ^ B vanish due to conflicting (5-function requirements. 
As usual we can assume N at rest with respect to the plasma, P = (M, 0,0,0), and defining 
c = cos 9 as the angle between pl and k we have 

A a (P L + K )= ^ = M(Kg '_ lKM , A (P H - K) = M(K -^ KM . (59) 

The loop integral of the first thermal term can be reduced to an integral over c and z = 
2\k\/M > as follows: 

f 1 dc 

x (1 _ c2)2 - d/2 if(Mz/2, k) + f(-Mz/2, k)). 

The last two terms in (58) can be similarly computed, by shifting the integration variable K 
in order to obtain the same factor S(K 2 )n B ^(K ) for all three factors. By dropping the T = 
contribution in jV \ we obtain 

r" fM) „ r o a 2 {e-^M 2 /?f 2 - 2 r d r x dc (zm d - 3 x 

1 vertex \ M ) ~ 1 tree 1c _a n/'J In 1\ / aZ flC fi „2VW/2 X 



167T 2 r(tZ/2-l) J (l- C 2)3-rf/2 

32 . . , 1 + c 2 
-^n B (Mz/2)-8— 2 



x { ^n B {Mz/2) - 8— -{n B {Mz/2) - n F {Mz/2)) \ . (61) 
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This expression has collinear divergences at c = ±1 (vector collinear with H or L) and soft 
divergences at z = (vector with vanishing energy). There are no divergences at z = 1; the 
apparent ones are treated taking the principal part. The structure of the singularities is plotted 
in fig. 4b. 

Turning now to the self-energy contributions, we observe that they can be written as in (29) 
with mjj top replaced with m 2 H jgauge + TO i, gaU ge an d — 1 with Zf^ ne ° + Z F — 2. The definition 
of the H thermal mass and wavefunction renormalization given in section 3.3.1 also holds in the 
gauge case (of course substituting lit with the gauge contribution to the H self-energy). We 
obtain Zf^ ugc = 1 and recover the known value of the thermal mass, rn 2 HgSMgc = g 2 T 2 /A. The 
calculation of the L self-energy requires more care. Like in the top case, what actually enters 
in the calculation is not the full L self-energy, but only its real part, which we call S(P^): 

E(P L ) = i 9 y- d J 0^A B (K)A F (P L + Kh,(If L + t)l, 

9^,4-^/9 W [ d d K {$ + If L )n B {K) + $n F {K) 
= 2g p (d/2 -I) J j^y d 2n5(K ) Jj^K) 2 

= g 2 ^- d (d/2-l) J j0zTm 2 )j^[nB(K o )+n F (K o )}, at P t 



0. 



The thermal mass can be defined again as the pole of the corrected propagator, \jyP ' F — S): 
squaring [f L — S) 2 ~ P 2 — {j^, S} we obtain the same expression as in eq. (18b): 



jgauge = \r L ^} = 2g 2 J 0^5(K 2 )[n B (K o )+n F (K o )} = ^T 2 . (62) 



Zl is the residue at the pole, explicitly given by Z F = 1 — Tr(7 S)/4Pi in the rest frame of 
the plasma, such that (see e.g. [13]): 

4o 2 / M\ d//2 ~ 2 dll — 1 r°° 
Zl = 1 ~ lt~ 2 (j J dj2~2 J, dZ zd ' 3 MMz/2) + n F (Mz/2)). (63) 

Z F contains infra-red divergences, which emerge when the loop momentum k is parallel to p F 
and are therefore of collinear type. 

3.4.2 Real gauge thermal corrections 

As already mentioned, the real corrections (diagrams b3, c3, d3 and d4 in Fig. 6) take the 
form (26). To simplify the phase space we can proceed like in the top case and repeat the steps 
in Eqs. (36), (37), (40), (39) and (41) with the substitutions P Q -> K and P v -> P H . In the 
gauge case we obtain 

i yy |2 9X 2 2 2(1-^ + ^(1-^/2-1) 

l^gaugel = 2A g z*y(l-y) " ( } 

Notice that |^g a u g e| 2 has soft (z = 0) and collinear (y = ±1) singularities. These lead to IR 
divergences in the real corrections as shown in fig. 4b. Like in the top case, we now express the 
real corrections as the sum of the following four integrals over y and z, corresponding to the 
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four regions in fig. 4: 



Vn^lha = rl ee ^^dy^dz(l-n F (P L0 ))(l+n B (K Q ))(l + n B {P^ 

2 roc poo 

T NH ^ LA = r° tiee ^j^ dyj^ dz(l-n F (P L0 ))(l+n B (K ))n B (P° H )F g (y,z), (65b) 

TnA^LH = r tree^ / / dz{l - U F (P L0 ) )n B (K ) (1 + U B (P°) )F g { V , z) , (65c) 

2 /-0 /-oo 

r^^A = -r t ° rccT f^ / / ^MPlo)(i + M^o))(i + ^(p£))p 9 (^), (65d) 

1D7T J_ OQ Ji 

where 

M r M n M , 

= y[l + *(y-l)], ^o = y^, P£ = y (l-^), (66) 

p , , = „ (e-^/jif- 2 |y(l - y)(l - ^)| d/2 ~ 2 [2(1 - g) + (1 - y)z\d/2 - 1)] 
j " r(d/2-l) y{l-y)z^\z\ i - d 

Real corrections harbour collinear divergences from y = and y = 1 and soft divergences from 
z = 0. Apart from the usual infra-red divergences (arising from soft and collinear configurations 
where |^g aU ge| 2 is singular), there are new purely thermal infra-red divergences, arising from 
configurations where a boson has vanishing energy such that their bosonic thermal distribution 
n B diverge, as anticipated at the beginning of section 3. 

3.4.3 Sum of the virtual and real contributions 

In principle, we should now proceed to compute explicitly the integrals in (61), (63) and (65). 
In practice they are so difficult that also checking the cancellation of IR divergences would need 
ad hoc tricks. In order to proceed in a systematic way, from now on we focus on the relevant 
limit T <C M, dropping all exponentially suppressed terms. 

For the real corrections we use the same observation already discussed in the top case: the 
dominant terms suppressed by powers of T/M arise from configurations where particles have 
energy ~ T <C M; and kinematics of 2 -H- 2 scatterings and 3^1 decays allows only one 
particle in turn to have small energy, such that we only need terms with one power of the 
thermal distributions in (65). In this way we split the thermal real contributions in three parts, 
K ~ 0, Pfj ~ and P£ ~ 0. 

It turns out that the divergences of the real contributions from the Kq ~ region cancel 
those of the vertex contribution, eq. (61), and that the divergences of the real contributions 
from the Pjy ~ and P° ~ regions cancel divergences from Z^. 

We therefore divide the calculation in these two groups to show the cancellation of diver- 
gences as soon as possible and to handle simpler expressions. 

• Kq ~ region and vertex contributions. The real contribution from the Kq ~ region 
emerges from (65a) and (65c) and corresponds to z ~ 0. Here we follow the same 
technique as in the top case: we expand F g (y, z) in powers of z including only those few 
powers that can modify the result up to 0(T 2 /M 2 ). We find 

9 2 f^ /E/2 M\ d -*T(d/2-2)T(d-4) [ 



- 



ro V_ „-7i5/2- 

tree 87r 2 V p. J r(d-3) 



2T 

M 



+ ( j^J (d-4)(d-3y(d-2)U d ^ 2 (l) 



Lid_ 4 (l)+ 

(68) 
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where Li n (z) is the polylogarithm. Expanding in d = 4 — 2e around d = 4 one gets a 
variety of infra-red divergent terms of the form 1/e 2 , 1/e, (lnT)/e, T 2 /e, InT, In 2 T as well 
as finite terms of the form T 2 and T 2 In T. 

A similar approach allows to evaluate the vertex contribution. The real Kq ~ and vertex 
contributions are separately lengthy expressions containing soft and collinear divergences, 
which manifest as poles at d — 4; however, a spectacular cancellation takes place and their 
sum is finite and simple 4 : 

n 2 + In 2 9 T 2 

r^— 3— a 2 W2 - (69) 

• Pjj ~ and Pi ~ regions and contribution. The real contributions from the P% ~ 
and P° ~ regions emerge from (65a)-(65b) and (65a)-(65d) respectively and correspond 
to the points (y, z) = (1, 1) and (y, z) = (0, 1) respectively. Thus, as in the top case we can 
now perform a Taylor expansion around those points. The result is rather complicated 
and contains a pole at d = 4 but its sum with the Z L contribution gives again a simple 
finite result: 

n 1/2 + In 2 0T 2 

Summing now (69) and (70) with the thermal mass contribution we obtain the result in (55): 
n 2 + ln2 9 T 2 n l/2 + ln2 9 T 2 n 1 9 T 2 

tree 3 «/ M 2 tree 3 «/ M 2 tree 2 y M 2 U ' ^ ^ 

4 Conclusions 

We performed a full NLO computation of one of two key ingredients for thermal leptogenesis: 
the interaction rate of the right-handed neutrino. This unifies partial contributions already 
included in previous works (2 — > 2 scatterings, thermal masses) and gives a very simple result, 
that we can summarize in one line: 



T N (E = M) 



M\ 2 (M) 



57T 



29 . . 21A 2 , T 2 

1 + — (3a 2 + « y )-^-^^ 

V * " 0.5(0) 

3%(2.5%) 5%(2%) 



(72) 



up to terms suppressed by higher powers of T/M or higher orders. At the same order, its 
thermal average 7^ is given by eq. (25). This is the only interaction rate that enters in the 
Boltzmann equations for the right-handed neutrino abundance and includes all the various de- 
cay and scattering terms considered in previous papers. A(M) is the neutrino Yukawa coupling 
renormalized at the right-handed neutrino mass M. 

The numerical values renormalized at the weak scale (at 10 10 GeV) are shown, having fixed 
Aft = (rrih/v) 2 assuming = 125 GeV. The total correction to the final baryon asymmetry 
ranges between 0.5% at 10 10 GeV and —3.5% at the weak scale, slightly depending on the value 
of rhi. 

The main difference between our result and previous results is the absence of infra-red 
enhanced corrections of order g 2 In rrih/M (where ~ gT is the thermal higgs mass), which 

furthermore eq. (68) has a pole at d = 5 because of the singular behavior of Li ra (l) for n — > 1. This emerges 
at finite temperature when the purely thermal singularity, Ub{E) for E — > 0, hits those that are also present at 
T = 0. This pole, as well as the poles at d — 4, cancels when adding up with the vertex contribution. 



19 



cancel out as dictated by the KLN theorem after performing our full quantum and thermal com- 
putation at NLO. In practice this means that 'gauge scatterings' and 'higgs scatterings' must 
be removed from codes for leptogenesis [5]: this makes computations simpler and effectively 
reduces 77V, thus enhancing the efficiency of leptogenesis 77 oc 1/jn- 

Furthermore, 7^ crucially enters also in the Boltzmann equation for the lepton asymmetry. 
In order to get the whole leptogenesis at NLO a second computation is needed: NLO corrections 
to the rate of CP-violating processes, which have been addressed at leading order in e.g. [14, 5, 
15]. This second more complicated computation is not performed in this paper; presumably it 
can be performed with the same techniques and the result has the same form, with cancellation 
of infra-red divergences. If the quantum version of Boltzmann equations will turn out to be 
needed [16] (which could be possibly the case in the presence of flavor oscillations), our precise 
computation of 7^ will remain needed. 
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Erratum 



Our computation in [1] has been confirmed by [2] up to a difference in one equation: when computing the 
reduction in the N — > LH rate due to thermal masses of L and H we incorrectly just considered the reduction 
of its phase space, neglecting that its amplitude is not constant. The full correct formula had been given in 
eq. (89) of [3]: 

T N (E-M) * —(1-2 — ). 

This means that below eq. (18) of [1] K T = -(2m| + m 2 H )/T 2 gets replaced by K T = -2m 2 H /T 2 . After this 
correction equations (27) and (28) read: 

^higgs = _j K to P = Q 
and the coefficient of the thermal effect T 2 /M 2 in the final result in equation (72) becomes 



T N (E = M) 



M\ 2 {M) 



29 , n . 21Af , T 2 

1 + 32^ (3a2 + av) - 32^ " Xh AP 



8tt 

in agreement with [2], These corrections have been fully implemented in the previous text [4]. 
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